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ABSTRACT 

We present the first observations of propane (CaHg) on Titan that unambiguously resolve propane 
features from other numerous stratospheric emissions. This is accomplished using a. R — Xj 5\ ^ 10^ 
spectrometer (TEXES) to observe propane's 1^26 rotation-vibration band near 748 cm^^. We find a 
best-fit fractional abundance of propane in Titan's stratosphere of (6.2 ±1.2) x 10^^ in the altitude 
range to which we are sensitive (90-250 km or 13-0.24 mbar). 

Subject headings: planets and satellites: Titan, infrared: solar system, molecular data 



1. INTRODUCTION 

Titan's thick atmosphere is simultaneously analogous 
to and extraordinarily different from that of Earth. Both 
atmospheres are composed mostly of nitrogen. Both at- 
mospheres have a major component that can exist as a 
solid, liquid, or gas (water on Earth, methane on Titan). 
Both have similar vertical temperature-pressure struc- 
tures, although Titan's atmosphere is ~200 K cooler and 
is greatly extended relative to Earth due to much lower 
surface gravity. On Titan a complicated network of pho- 
tochemical reactions leads from methane (CH4) to the 
formation of numerous heavier hydrocarbons, including 
propane (CsHg). Measuring the abundances of Titan's 
hydrocarbons provides an important test of models of Ti- 
tan's atmospheric chemistry and, more generally, of our 
fundamental understanding of atmospheric chemical and 
physical processes. 

Photochemical models (e.g. lYung. Allen^ fc Pintd 
rm84t ITo ubla,nc et al.1 1199^ IT.ara et abj FmM 
iBanaszkiewicz et al.l l2000j) predict that propane is 
formed via C2H5 + CH3 M ^ CsHg + M where CH3 
is created by methane photolysis and C2H5 by C2H4 -f 
H + M — * C2H5 -I- M. The primary removal mechanism 
for propane is condensation at the cold tropopause, 
although photolytic destruction and reaction with C2H 
co nsume ^5-1 5 % of the propane formed. In the model 
of ILara et all l)1996|) propane condensed onto Titan's 
surface accounts for ~1% of the carbon removed from 
Titan's atmosphere. 

Detections of propane's z^26 rotation-vibration band 
near 748 cm~^ were based on Voyager I and Infrared 
Space Observatory (ISO) spectra. Voyager spectra had 
4.3 cm-i resolution {R = X/5\ = 170 at 748 cm"' 
and initial claims of a propane abundance of 2x10 
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Maguire e t al."1981') were later reduced to 7±4 x 10~^ 
Cou stenis et al. 198^ . The resolution of these spec- 
tra left every band of propane blended with other emis- 
sion features. At 748 cm~^ the V2q band of propane 
was blended with features of acetylene (C2H2) and the 
derived abundance (or even detection) of propane is 
strongly dependent on the assumed C2H2 vertical abun- 
dance profile. Spectra taken by ISO were of higher res- 
olution (i? « 1850 at 748 cm~^) but were still unable 
to resolve propane fr om acetylene. From these data 
ICoustenis et al.l ()2003j) derived a nominal propane abun- 
dance of 2±1 X 10~^, although they acknowledged the 
substantial difficulties of disentangling propane emission 
from that of acetylene. 

2. OBSERVATIONS 

Observations were made at the 3.0 meter NASA 
Infrared Telescope Facility (IRTF) using the TEXES 
mid-infrared high-resolution grating spectrograph 
on 13, 14, and 20 December 2002 



l|La,cv et a,1.l l2f 

UT in several overlapping spectral settings around 
propane's fundamental vibration at 748 cm^^. 

At this wavenumber a single spectral setting covers 
~5 cm^^ at a resolution of R k, 10^ with slight gaps 
between spectral orders. The strongest ^26 propane 
lines are blocked by telluric absorption at 748.3 cm~^. 
The diffraction limit of IRTF at these wavelengths and 
Titan's diameter are both w 0'.'9. With a slit width 
of 1'.'5 our data represent a disk average spectrum of 
Titan. The reduced spectrum from each night, corrected 
for the appropriate Doppler shift, is shown in Fig. ^ 
along with the final combined spectrum. All data 
were reduc ed using the stan dard pipeline reduction 
detailed in ILacv et all 1)20021) . Telluric transmission 
was corrected with Callisto on t he first two nights and 
theBecklin-Neugebauer object l)Becklin fc NeugebaueJ 
Il967^ on the third nig ht. We flux cali brated using 
observations of (3 Gem ()Cohen et alJ '1995l on the final 
night, and estimate that this calibration is accurate to 
10-20%. 

3. RADIATIVE TRANSFER MODELING 

We use a new line-by-line code (|Rod l2002|) , divid- 
ing Titan's atmosphere below 1000 km into 50 spherical 
shells evenly spaced in log (P) and using spectral bins 
of 3 X 10~^ cm~^ in order to resolve even the narrow- 
est emission lines. Line parameters for acetylene (C2H2) 
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Fig. 1. — Reduced spectra from each night (offset for clarity) and 
a final combined spectrum are shown on the upper panels (a,c), 
while the calculated telluric transmission ( Ro#2003, Appendix A) 
is shown on the lower panels (b,d). When compared with an empir- 
ical absorption coefficient spectrum of propane (shown at bottom 
of a,c in arbitrary units; See § 3 for more details.), nearly all un- 
labeled features in the observed spectra are seen to be associated 
with propane. 



and HCN are from HITRAN l|Rothman et al. 1998|) . 
with the HCN hne positions adjusted to ag ree with 
the observations of iDuxburv fc Gan j l|1989() . The 
temperature-pressure profile is the 'recommended' pro- 
file of Yclle et al. (1997). Scattering effects are ignored 
and Titan's haze is modeled as a single layer with a lower 
cutoff of 50 km and with haze opacity scale height equal 
to the gas density scale height. The abundances of HCN 
and C2H2 are each parameterized with their fractional 
abundances at 1 mbar {FHCN,im,bar, Fc2H2,imbar) and 
with the slope of /o(/(abundance)-/o(7 (pressure) {riHCN, 
nc^H-i), i.e. Fhcn = -FffC7V,imfcar(^'mfcar)""^"- Wc in- 
vestigated both a constant propane vertical profile as 
well as scaled versions of an abundance profile pre- 
dicted by photoche mical modeling (Laraetal. 1996; 
iBanaszkiewicz et al.l 120011) . Species are held to satura- 
tion vapor pressures below their condensation altitudes. 
In order to calculate Titan's total fiux, the model is run 
at 32 points from the center of Titan's disk to off the 
edge of Titan's solid limb, with the calculation points 
more closely spaced near the limb. 

The modeling of propane's emission spectrum re- 
quires additional discussion. Several linelists exist for 
propane's V2& ba nd at 748 cm~^. These inclu de: the 
GEISA databank ij.Tacauinet-Husson et al. 1999() : an un- 
published list based on fits to t he laboratory s pectra 
of iNadler fc Jennings! 1)19891) and iHillman et aD (|1992D 
(Blass et al. 1988); and a fist generated using the code 
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Fig. 2. — (a ) Comparison betw een low resolution laboratory 
measurements of lGiver et al.1 ITflSl^ at 200 K with calculated spec- 
tra for the same conditions using the available line lists. Calculated 
spectra are scaled to match the observed spectrum as closely as 
possible and offset for clarity. There is good agreement in the 
wings of the band, however over ~740-760 cm~^ there is very poor 
agreement between calculations and observation, (b) Comparison 
of calculated spectra using the available line lists at 175 K and 3.2 
Torr with an observed transmission spectrum of propane at the 
same conditions from Hiilman et al. 1 1992) . There is good agree- 
ment at some wavenumbers (e.g. 748.03-748.15 cm~^) and poor 
agreement at others (e.g. < 748.0 cm"-*^). 



of iTvpkd l)1976f) and the parameters of iGassler et alJ 
l)1989() . Figure El shows calculated spectra using each 
of these li nelists com pared to low-resolution laboratory 
spectra of iCiver e^ al. ( 1984) a nd to the high-resolution 
spectra used in iHUlman et all l)1992|) . Given the poor 



fit of the calculated spectra to the observed spectra in 
Fig. [3 at the wavenumbers we observed, we opted to use 
the low-pressure (1-3 Torr) low-temperature (~175 K) 
high-resolution (Kk, 3 x 10^) transmission spectra of 
IHillman et aD l(l99^ ^ to construct an empirical absorp- 
tion coefficient [kv) spectrum for propane. This involved 
fitting for the smoothly varying baseline and then con- 
verting to fcjy units (cm^^ amagat"^) using the sample's 
known pathlength (30 cm), pressure, and temperature. 
In using this empirical spectrum we are implicitly 
assuming that k^ is independent of temperature. Our 
propane observations primarily probe regions of Titan's 
atmosphere that are at 135 to 175 K and are therefore 
reasonably well matched by the conditions of these lab- 
oratory data used to construct k^. Most of our spectral 
fitting was performed using this empirical for propane. 
As a test of the importance of the mismatch between Ti- 
tan atmospheric and laboratory temperatures we re-fit 
our data using the Blass et al. linclist scaled to fit our 
empirical k^, at 1 75 K, along with the partiti on sum pa- 
rameterization of iFischer fc Gamachd l)2002|) . to repre- 
sent propane. 

4. FITTING MODEL TO DATA 
^ Available at |http: //diglib.nso.edu/nso-U^er.html| 



Titan's propane 
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Fig. 3. — Comparison of observed spectrum with model con- 
volved to the observed spectral resolving power. Spectra are offset 
for clarity and residuals are shown on the same scale. The best fit 
is obtained using our empirical ki, spectrum for propane. 



We fit model to data with an 'amoeba'-type mini- 
mization while varying six parameters: total haze opac- 
ity {Thaze), propane (FcsHg), acetylene iFc2H2Smbar, 
nc2H2), and HCN {FHCN,imbar, nncN)- Figure |3l com- 
pares the observed spectrum with the best-fit model with 
and without propane, in one case using the Blass et al. 
linelist and in the other using our empirical ki, to rep- 
resent propane. The detection of propane is unambigu- 
ous and the best-fit abundance assuming a constant ver- 
tical profile and using our empirical k^, is 6.2 x 10~^. 
For acetylene and HCN our best fits are characterized 
by FC2H2 = 4.4 X 10-6(P„fc,,)-0-45 and Fhcn = 3.8 x 
10-'^(P„6ar)"''-^*- The haze is best fit by Thaze = 0.21, 
although this value is extremely dependent on the as- 
sumed vertical haze profile and should not be over inter- 
preted. 

We found that reasonable fits to the HCN and C2H2 
features required a non-constant vertical abundance. 
Even with our two variable log-log parameterization we 
see non-zero residuals on these features, suggesting a 
more complicated vertical structure which we will explore 
further in future work with similar spectra covering many 
more transitions. Other significant non-zero residuals are 
explained as: an unidentified species at 747.35 cm~^, an 



imperfectly corrected telluric feature at 748.68 cm , not 
fitting for the i^CCHa/CzHa ratio at 748.88 cm^^, and 
an H^'^CN feature at 749.02 cm"-'^ that is missing from 
the model. 

Our focus here is propane and thus we do not present 
uncertainties on the vertical abundance parameters for 
HCN and C2H2. As is apparent in Fig. ^ random er- 
rors are negligible in our spectra. The other sources of 
uncertainty are: flux calibration, errors associated with 
our empirical absorption coefhcient method, the extent 
to which propane abundance is not constant with alti- 
tude, the extent to which Titan's te mperature-pressur e 
profile deviates from that derived bv lYelle et alJ ljl997fl . 
and the extent to which propane's abundance varies with 
latitude. 

The uncertainty in flux calibration is at worst 20%, 
which maps almost linearly into a 20% uncertainty in 
propane abundance. The primary potential error from 
our empirical absorption coefficient method is the mis- 
match of laboratory sample temperature (175 K) to the 
temperature of the propane in Titan's stratosphere. In 
our model spectra > 80% of the propane flux arises from 
the pressure range 13-0.24 mbar (90-250 km altitude), 
where temperatures range over 135-175 K. To investigate 
this uncertainty we re-fit using the Blass et al. linelist 
(See Fig. 131). The overall fit is lower quality than with 
our empirical k^, and the best-fit propane abundance is 
only l.libO.l times more. 

Recent photochemical m odels (e.g. ILara et aD Il996t 
IBanaszkiewicz et al.l I2000|) predict that the propane 
abundance is relatively constant over 400-850 km (11- 
1.5 /ibar). Above these altitudes propane is depleted 
by photolysis and reaction with C2H, while below this 
range the effect of eddy diffusion and the cold-trap of 
the tropopause is to reduce the propane abundance. Our 
observations are primarily sensitive to the pressure range 
13-0.24 mbar (90-250 km altitude), as shown in Fig.|^. 
We re-fit our spectra using the model predicted propane 
vertical profile and found a best fit when the model pre- 
dicted profile is multipled by a factor of 2.9. F igure Bb 
shows this best-fit scaled profile of ILara et a l. ( 1993) 
along with the best-fit constant vertical abundance pro- 
file. The quality of the fits to our data do not distinguish 
between these two cases.^ 

The uncertainty in the lYelle et alJ l)1997l) temperature- 
pressure profile is not well quantified, and therefore our 
reported uncertainties do not account for this source of 
error. Roughly, a change in stratospheric temperature of 
±5 K requires a =f20% change in propane abundance. 

There is no strong evidence that propane abundance is 
constant with latitude on Titan. A first attempt at mod- 
eling the season al variation of photochemistry on Titan 
ijLebonnois et al.-2001) predicts that propane abundance 
will vary by a factor of a few from north to south depend- 
ing on season. With TEXES on an 8-10 meter telescope 
we would be able to measure the abundance as a func- 
tion of latitude, which will place new strong constraints 
on models of seasonal photochemistry. 

5. SUMMARY 

We present the first spectrally resolved detection of 
propane in Titan's atmosphere. We measure propane's 
fractional abundance to be (6.2 ± 1.2) x lO"*", assum- 
ing the recommended temperature-pressure profile of 
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Propane fractional abundance 

Fig. 4. — (a) Contribution function of propane emission at the 
middle of the 748.10 cm~^ feature. The contribution function is 
the partial derivative of propane emission with respect to log P. 
(b) Comparison of the two best-fit vertical abundance profiles. 



The current theoretical understanding of propane's 
vibration-rotation spectrum is inadequate and avail- 
able linelists do not fit our Titan spectra nor labora- 
tory propane spectra well. Therefore we fit for Titan's 
propane abundance using an empirical absorption coeffi- 
cient spectrum derived from low temperature laboratory 
spectra of propane. 

The most important observational advance that can 
be made towards detecting new species in planetary at- 
mospheres is to increase the spectral resolving power to 
R >fewxlO'' such that individual lines are resolved and 
separated. High spectral resolution allows us to sepa- 
rate the contribution of a minor species (propane) from 
a strongly emitting species (acetylene). 



lYelle et aTl l|1997f) and a constant abundance with alti- 
tude and latitude. Our observations are primarily sen- 
sitive to 13-0.24 mbar (90-250 km altitude). Alterna- 
tively, our data r equire the predicted propane profile of 
iLara et aTl l|199fifl to be increased by a factor of 2.9. 
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